Role of electronic correlations in the Fermi surface formation of Na^CoC^ 
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Band structure of metallic sodium cobaltate Na^CCh (x=0.33, 0.48, 0.61 0.72) has been investi- 
gated by local density approximation+Hubbard U (LDA+C/) method and within Gutzwiller approx- 
imation for the Co-t2 g manifold. Correlation effects being taken into account results in suppression 
of the e' g hole pockets at the Fermi surface in agreement with recent angle-resolved photo-emission 
spectroscopy (ARPES) experiments. In the Gutzwiller approximation the bilayer splitting is signif- 
icantly reduced due to the correlation effects. The formation of high spin (HS) state in Co ci-shell 
was shown to be very improbable. 

PACS numbers: 74.70.-b; 31.15.Ar; 71.45.Gm; 71.10.-w 
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Introduction. Puzzling properties of sodium cobal- 
tate Na^C02 are the topic of many recent theoretical 
and experimental investigations. This material holds 
much promise for thermoelectronics due to its large ther- 
mopower— together with the relatively low resistivity 2 . 
The discovery of superconductivity with T c about 5K 
in Nao.33CoO2T.3H2O 2 revived the interest in lamellar 
sodium cobaltates. Moreover, the charge and magnetic 
long range orders on the frustrated triangular lattice of 
cobaltate is of the fundamental interest. The band theory 
predict the complicated Fermi surface (FS) with one large 
hole pocket around the T = (0, 0, 0) point and six small 
pockets near the K = (0,^,0) points of the hexagonal 
Brillouin zone at least for x < 0.5 ?,4; . However, intensive 
investigations by several ARPES groups reveal absence 
of six small pockets in both Na K C02-yH20 and in its 
parent compound Na^CO^^^. 

The disagreement between ARPES spectra and ab- 
initio calculated band structure points to the importance 
of the electronic correlations in these oxides. Other evi- 
dences for the correlated behavior come from the data on 
an anomalous Hall effect and a drop of the thermopower 
in holistic magnetic fielcU^. 

The six hole pockets are absent in the L(S)DA+£/ cal- 
culation o 12 ' 13 . However, in this approach, the insulating 
gap is formed by a splitting of the local single-electron 
states due to spin-polarization, resulting in a spin polar- 
ized Fermi surface with an area twice as large as that 
observed through ARPES. Moreover, the long range fer- 
romagnetic order has been set by hand because of lim- 
itation of LDA+f/. The predicted large local magnetic 
moments as well as the splitting of bands can be consid- 
ered as artifacts of the L(S)DA+[7 method. 

Although LDA+U method is usually applied to de- 
scribe insulators, there are some achievement in inves- 
tigation of metals and metallic compound a 14 ' 15 . To an- 
alyze the effect of electronic correlations on the Fermi 
surface formation in sodium cobaltate we employ both 
LDA+J7 method and a Gutzwiller approximation. 



Co d-level splits by crystal field of oxygen octahedron 
in lower ti g and higher e g bands. The deficiency of 
Na in Naa;C02 introduces additional holes in the sys- 
tem. Cobalt, having d 6 configuration and filled t2 g shell 
in parent NaC02, is nonmagnetic. But in nonstoichio- 
metric compound part of Co ions become magnetic with 
local moment about 1/j.b- This value is provided by d 5 
configuration and one hole in ti g shell. However, the ex- 
periments revealed the magnetic susceptibility at room 
temperature that is much higher than it was expected 
for dilute magnetic impurity in non-magnetic solvent. 
Explanation of this anomaly was suggested in Re£^ as 
transition from low-spin (LS) state with six ci-electrons 
on t2 g shell to high-spin (HS) state with five d-electrons 
on ti g shell and one electron on e g . The possibility of 
such transition will be discussed below. 



LDA+U results. Nao.6iCo02 crystallize in the 
hexagonal unit cell {PQ^/mmc space group) with 
a=2.83176(3)A and c=10.8431(2)A at 12 K Ref 18 . 
Displacement of Na atoms from their ideal sites 2d 
(1/3,2/3,3/4) on about 0.2A are observed in non- 
stoichiometric cobaltates for both room and low temper- 
atures. This is probably due to the repulsion of a ran- 
domly distributed Na atoms, locally violating hexagonal 
symmetry 18 . In the present investigation Na atoms are 
shifted back to their 2d ideal sites. In order to avoid 
charge disproportionation which can arises from some 
Na distribution if the supercell is used in calculation the 
change in the Na concentration has been considered in 
virtual crystal approximation (VCA) where each 2d site 
is occuped by virtual atoms with fractional number of 
valence electrons x and a core charge 11 — x instead of 
Na. Note, that all core states of virtual atom are left un- 
changed and corresponds to Na ones. We have chosen 4s, 
4p, and 3d states of cobalt, 2s, 2p, and 3d states of oxy- 
gen, and 3s, 3p, and 3c? states of Na as the valence states 
for TB-LMTO-ASA computation scheme. The radii of 
atomic spheres where 1.99 a.u. for Co, 1.61 a.u. for oxy- 
gen, and 2.68 a.u. for Na. Two classes of empty spheres 
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FIG. 1: (Color online) Band structure of Na^OoC^ for x equal 
to 0.33 (a), 0.48 (b), 0.61 (c), and 0.72 (d), obtained within 
LDA is shown by the black (dashed) curves. Band structures 
for the same doping concentrations within LDA+U are shown 
by the red (solid) curves. 



(pseudo-atoms without core states) were added to fill the 
unit cell volume. 

Crystal field of oxygen octahedron splits Co e?-band 
into 2-times degenerate e g and 3-times degenerate ti g 
subbands (without taking into account spin). LDA cal- 
culations shows that those manifolds are separated on 
about 2 eV— . Here partially filled ti g subband crosses 
the Fermi level whereas e g subband due to strongly hy- 



bridization with nearest oxygen atoms is positioned well 
above the Fermi level. The procedure proposed in Reft^ 
allows one to calculate Coulomb repulsion parameter U 
taking into account the screening of localized d-shell by 
itinerant s- and p-electrons. Resulting U is equal to 6.7 
cV. However, the presence of the ti g -e g splitting give 
the reason to take into account an additional screening 
channel provided by the less localized e g electrons. The 
value of U for ti g orbitals calculated using "constrained 
LDA" method 20 , in which the screening by e g electrons is 
also taken into account is equal to [7=2.67 eV. The later 
value was used in the present calculation for all doping 
concentrations x. Hund's exchange parameter J depends 
weakly on screening effects due to its "on-site" charac- 
ter. Its value was calculated within "constrained LDA" 
method and is equal to 1.07 eV. 

First, we have verified the possibility of HS state for- 
mation on Co d-shell. For this purpose the unit cell of 
Nao.6iC02 with two Co atoms was considered. We have 
started from a saturated A-type antiferromagnetic con- 
figuration with five electrons on the ti g and one on the e g 
shells. Small [7=2.67 eV does not stabilize such magnetic 
configuration and LS state was obtained. Increasing [7 
up to 5 eV however results in HS state with large local 
magnetic moment about 1.96 \ib- Nevertheless, this HS 
state has the total energy about 1.75 eV higher then the 
energy of a LS state. This large difference in total en- 
ergy of both considered spin states arises form hexagonal 
structure of cobaltates where the angle of Co-O-Co bond 
is close to 90° in contrast to almost 180° for, e.g., RC0O3 
(R=La, Ho) . In the later case the e g band has the width 
about 3-5 eV and its bottom lies just above the Fermi 
level. The system wins energy of 2 J forming a HS state 
overcoming the gap energy which is less than 1 eV. Due 
to this fact difference between LS and intermediate spin 
states in RC0O3 is less then 250 meV— . The angle of 
Co-O-Co bond is close to 90° in cobaltates and results in 
a weak overlap between e g orbitals and hence in a narrow 
e g band with larger gap between it and the ti g band. Our 
calculation confirms that formation of HS in Na^CoC^ is 
rather improbable and cannot be stabilized by any dis- 
tortion of crystal structure or clusterization proposed in 
Reffii. Local magnetic moments on Co sites can arise 
only due to holes doping due to Na atoms deficiency. 
Those holes order on Co atoms and form nonmagnetic 
Co 3+ and magnetic Co 4+ ions with d 6 and d 5 configura- 
tions, respectively. In the further calculations only the 
LS state was considered. 

The ordering of holes on <2g shell and correspond- 
ing long-range magnetic and charge order for Nao.5Co02 
arise probably due to specific arrangement of Na atoms. 
This arrangements were observed experimentally 2 ^ for 
several doping concentrations including x = 0.5. Proper 
description of such order within "unrestricted Hartree- 
Fock" gives strong spin and orbital polarization and lo- 
cal magnetic moment of about l/i^ on Co 4+ sites as 
well as insulating ground state with sizable gap. To de- 
scribe non-ordered systems the implementation of "re- 
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stricted Hartree-Fock" method is more suitable. In the 
later starting from non-magnetic configuration of d-shell 
with equal number of spin-up and spin-down electrons 
LDA+t/ method gives non-magnetic solution without 
spin or orbital polarization. Note, that the gap does not 
open and Na^OoC^ remains metallic for all Na concen- 
tration. 

Obtained band structure of Na^OoO^ for x=0.33, 0.48, 
0.61, and 0.72 are shown in Fig. [TJ Dashed (black) lines 
correspond to LDA results whereas solid (red) lines are 
the bands obtained by LDA+J7 method. Cobalt d and 
oxygen p states are separated by a small gap of about 
-1.25 eV for x—0.61 and £=0.72. However, this gap dis- 
appears for lower doping concentration since the d band 
goes down when the number of d electrons decreases. 
The presence of the two C0O2 layers within a unit cell 
due to alternation of the oxygen arrangement results in 
a bonding-antibonding (bilayer) splitting, also present in 

Fig. in 

The degeneracy of the ti g levels is partially lifted by 
the trigonal crystal field distortion which splits the for- 
mer into the higher lying a\ g singlet and the lower two 
e' g states. However, slight difference in occupation num- 
bers of ai g and e' g orbitals (0.714 and 0.886 respectively 
for .x=0.33) results in significant difference in LDA+C7 
band structure in respect to LDA one. The energy of the 
less occupied a\ g orbital increases for both spins, whereas 
all e' g bands go down (the total a\ g -e' g splitting becomes 
0.21 eV for ir=0.33). This makes six e' g Fermi surface 
hole pockets to disappear for small x values. Note, for 
all Na concentration LDA+J7 predicts large a± g pocket 
of Fermi surface centered around the T point in excellent 
agreement with ARPES spectra for x < 0.7. The ad- 
ditional electron pocket close to the T point appears in 
both LDA and LDA+J7 methods for a large doping con- 
centrations. It was discussed in our previous wor k 23 ' 24 
in connection with the electronic theory for the itinerant 
magnetism of this compound. 

Gutzwiller approximation. The free electron Hamilto- 
nian for Co02-plane in a hole representation is given by: 
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FIG. 2: (Color online) Band structure of Na^CoCh for x 
equal to 0.33 (a), 0.48 (b), 0.61 (c), and 0.72 (d), obtained 
in LDA+U is shown by the dashed (red) curves. Dispersion 
within the Gutzwiller approximation is shown by the solid 
(black) curves. 
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where df aa (df aa ) is the annihilation (creation) operator 
for the ti g hole at Co site f , spin a and orbital index a, 
nfaa = dl acr df acr , and t^ is the hopping matrix element, 
e a is the single-electron energies, and n is the chemi- 
cal potential. All in-plane hoppings and single-electron 
energies were derived previousl y 23 ' 24 from the ab-initio 
LDA calculations. Since LDA-calculated hoppings and 
single-electron energies do not depend much on doping 
concentration, we used here parameters for Nao.33Co02 
form Table I of Refi^i. To take bilayer splitting into ac- 
count, we also consider here the hoppings t"' 3 between 
adjacent Co0 2 planes. Their values were also derived 
from LDA results and are equal to t a z l9<lls = -0.0121, 

t^ 16 ' 91 = 0.0080, and = -0.0086. 



The Gutzwiller approximatio n 25 ' 26 ' 27 for the Hubbard 
model provides a good description for the correlated 
metallic system. Its multiband generalization was for- 
mulated in Ref.— . In this approach, the Hamiltonian 
describing the interacting system far from the metal- 
insulator transition for U 3> W . J = 0: 



H — H + E U a rif a ^nf a i, 
f ,a 



(2) 



is replaced by the effective non-interacting Hamiltonian: 

H -ff = - E ( £Q + S£aa - /*) nf «- 



f,a,<7 



E E^4*a^+c. (3) 
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Here, = tf£ ^Jqaa^qpa is the renormalized hopping 
between two lattice sites connected by the spatial vector 

(f - g)> Qa<r = ]_-n aa , n aa = (* | nfaa |*o) = (™faff) is 

the orbital's filling factors, x — 1 — ^2 n aa is the equation 

for the chemical potential. 5e aa are the Lagrange multi- 
pliers yielding the correlation induced shifts of the single- 
electron energies. The constant C — ^£ a<Tn aa was 

determined from the condition that the ground state en- 
ergy is the same for both Hamiltonians, ("Pol Heff l^o) = 
(tyg\ H \^f g ), where |^o) is the wave function of the free 
electron system ((3]), and \^ g ) is the Gutzwiller wave func- 
tion for the Hamiltonian @. The Lagrange multipliers 
are determined by minimizing the energy (&o\ H e ff \^o) 
with respect to the orbital filling factors n aa . This re- 
sults in the following expression for the single-electron 
energies renormalization: 

^ = ^E^(04. (4) 

It is this energy shift that forces the e' g bands to sink 
below the Fermi energy^. This is clearly seen in the 
doping-dependent evolution of the quasiparticle disper- 
sion within Gutzwiller approximation in Fig. O To ob- 
tain these figures we solved self-consistently Eq. (g]) and 
the equation on the chemical potential p.. 

The comparison of the Gutzwiller approximation re- 
sults with the LDA+C/ dispersion reveals few very inter- 
esting conclusions. First, both approximations result in 
a suppression of e' g hole pockets of the FS. Second, the 
bilayer splitting is strongly doping dependent and signif- 
icantly reduced for Gutzwiller quasiparticles in compar- 
ison with the LDA+U quasiparticles because the renor- 
malization coefficient, ^qaa^/Qfla, occurs not only for the 
in-plane hoppings, but also for the interlayer hoppings 
t^f . Third, when both bonding and antibonding ti g 
bands do not cross the Fermi level around the T point, the 
FS crossings are the same in both approximations (see 
Fig. [5^). It is the simple consequence of the Luttinger 
theorem which holds for both approaches. But for large 
x due to the larger bilayer splitting in LDA+/7 the Fermi 



surfaces become different, while the Luttinger theorem is 
still preserved. With increase of the doping concentration 
x the bandwidth of the Gutzwiller quasiparticles becomes 
closer to the LDA+C/ ones because band renormalization 
factor ^/q a <r^q_i3a comes closer to unity. 

Dynamical Mean Field Theory (DMFT) 
calculation o 31 ' 32 show that for the small U and non-zero 
J, e' g FS pockets can by stabilized. On the other hand, 
DMFT calculations of ReL^S confirms results of the 
Gutzwiller approximation provided that U and the 
crystal filed slitting are large. 

Conclusion. In the present work employing ab-initio 
"constrained LDA" method we obtained Coulomb repul- 
sion parameter [7=2.67 eV for t 2g orbitals taking into 
account the screening by e ff -electrons in addition to the 
screening by itinerant s- and p-electrons. Hund's exch- 
nage parameter was found equal to J=1.07 eV. 

Also, we have shown that due to the Co-O-Co bond 
angle close to 90° in Na2;Co02 the energy gap between 
LS and HS states is too large to be overcome by clus- 
terization or reasonable distortions of crystal structure. 
Thus we conclude that realization of HS state is highly 
improbable in these particular substance. 

Then, to analyze the effect of electronic correlations on 
the Fermi surface topology of NazCo02 we use two ap- 
proaches, non-magnetic LDA+[7 and Gutzwiller approx- 
imation for the Hubbard-type model based on the LDA 
band structure. Within LDA+C7 the energy of the less 
occupied a\ g orbital increases for both spins, whereas all 
e' g bands go down. This makes six e' g FS hole pockets to 
disappear for small x values, in agreement with ARPES 
for x < 0.7. Gutzwiller approximation also resulted in a 
suppression of e' g hole pockets at the FS. Moreover, the 
bilayer splitting was found to be strongly doping depen- 
dent and significantly reduced for Gutzwiller quasiparti- 
cles in comparison with LDA+[7 quasiparticles. 
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